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ABSTRACT 

Swami,  Shanmugam  Arumuga.  Ph.D.,  Purdue  University,  January  1969. 
The  Response  of  Bituminous  Mixtures  to  Dynamic  and  Static  Loads  Using 
Transfer  Functions .  Major  Professor:  William  H.  Goetz. 

The  applicability  of  the  concept  of  transfer  functions  for  bituminous 
concrete  has  been  examined  in  this  study  by  treating  the  test  specimen  as 
a  dynamic  system  in  which  the  input  is  the  applied  force  and  the  output 
is  the  resulting  displacement.  The  transfer  function  was  derived  from  the 
frequency  spectrum  of  the  material  using  a  simple  curve  fitting  technique, 
and  the  spectrum  itself  was  obtained  from  a  series  of  sinusoidal  load  tests 
of  varying  frequencies. 

The  transfer  function  derived  for  the  asphaltic  concrete  was  found 
to  be  of  form 


A(s  +  a1)(s  +  a2)(s  +  a^ 
G(s)  =  (s  +  b1)(s  +  b2)(s  +  b3)(s  +  bk) 


where  s  is  a  complex  variable,  A  is  a  constant  and  a  's  and  b.  *s  are  the 
roots  of  the  numerator  and  denominator,  respectively.  These  roots  were 
obtained  directly  as  corner  frequencies  in  the  asymptotic  approximation 
of  the  frequency  spectrum. 

The  scope  of  the  study  included  dynamic  and  static  testing  of  speci- 
mens of  asphaltic  concrete  of  two  different  compositions,  cut  along  three 
mutually  perpendicular  directions  from  laboratory  compacted  beams  with 
all  testing  being  carried  out  at  three  different  temperatures. 


xi 

The  follov/ing  primary  conclusions  may  be  stated  as  a  result  of  this 
study : 

1.  The  viscoelastic  or  time-dependent  characteristics  of  a  bitumi- 
nous concrete  can  be  represented  by  its  transfer  function  which  is  unique 
for  the  material  at  a  given  temperature. 

2.  The  roots  of  the  denominator  of  the  transfer  function  of  the  as- 
phaltic  concrete  tested  were  observed  to  be  real  and  distinct  which  in- 
dicates that  bituminous  concrete  behaves  as  an  overdamped  system. 

3.  Temperature  is  the  one  single  factor  which  most  affects  the 
transfer  function  for  an  asphaltic  concrete  material.  Increase  in  tem- 
perature increases  the  value  of  the  constant  A  in  the  transfer  function 
equation. 

U.     The  transfer  function  is  observed  to  be  a  powerful  tool  in  pre- 
dicting the  displacement  under  an  applied  load,  dynamic  or  static.  By 
treating  the  static  load  as  a  step  function  of  time,  the  static  load  dis- 
placement was  calculated  by  means  of  the  transfer  function  derived  ex- 
perimentally from  the  dynamic  test,  and  compared  with  the  measured  static 
load  displacement  value.  Excellent  agreement  was  noted  between  the  calcu- 
lated and  measured  values  of  the  displacements.  This  validates  the  con- 
cept that  the  transfer  function  represents  a  material  property  which  is 
independent  of  the  type  of  input . 

5.  By  the  use  of  the  transfer  function  and  without  assuming  any 
spring-dashpot  model,  it  is  possible  to  represent  the  time-dependent  be- 
havior of  asphaltic  concrete  by  a  fourth  order  linear  differential  equa- 
tion with  constant  coefficients.  The  coefficients  can  be  computed  from 
the  roots  of  the  denominator  and  numerator  of  the  transfer  function. 


INTRODUCTION 

It  is  only  in  recent  years  that  the  role  of  bituminous  concrete  as  a 
load -distributing  layer  in  flexible  pavements  has  been  well  recognized, 
even  though  it  has  been  used  in  the  upper  layers  of  the  pavement  for  a 
long  time.  The  results  of  the  AASHO  Road  test  (l)  demonstrated  that  re- 
duction in  pavement  deflections  and  improved  performance  were  effected 
through  the  use  of  asphalt -treated  materials  in  the  load-distributing 
layers.  In  these  applications  both  the  bituminous  material  and  the  aggre- 
gate serve  as  structural  elements . 

Currently  used  mix  design  methods  (2)  employ  stability  numbers  that 
may  provide  some  measure  of  quality  control  but  do  not  relate  mixture  con- 
stituents to  loads  imposed  on  the  pavement.  Further,  these  methods  do  not 
measure  any  fundamental  parameters  which  can  characterize  the  mechanical 
behavior  of  the  mixture  under  a  load,  static  or  dynamic,  nor  do  they  con- 
sider the  time -dependent  characteristics  of  the  bituminous  mixture. 

Placed  in  the  top  layer  of  a  pavement,  the  bituminous  mix  derives  its 
support  from  the  underlying  base  and  subgrade  layers  and  together  they 
form  a  composite  pavement  structure.  Hence,  an  essential  first  step  in 
the  rational  design  of  a  pavement  system  is  the  ability  to  estimate  the 
stresses  and  strains  induced  in  its  different  layers  under  a  given  type 
of  loading,  at  any  given  time,  through  a  knowledge  of  the  material  proper- 
ties determined  from  simple  laboratory  tests  that  will  have  significance 


Figures  in  parentheses  refer  to  references  in  the  Bibliography. 


under  unknown  and  very  complicated  field  conditions.  A  necessary  second 
requirement  is  the  availability  of  some  failure  criteria,  in  terms  of 
either  stresses  or  strains  for  the  material  of  each  layer,  so  that  the  de- 
signer may  fix  allowable  values  of  these  quantities.  It  is  obvious,  how- 
ever, that  in  the  absence  of  a  theory  or  a  method  which  gives  the  stresses 
and  strains  in  a  pavement  structure,  the  second  requirement  cannot  be  of 
much  practical  value. 

Use  of  Elastic  Theory 

Elastic  theory  has  long  been  used  to  investigate  the  stresses  and 
displacements  in  layered  systems.  Solving  for  the  stresses  in  a  pavement, 
Westergaard  (3)  assumed  the  pavement  to  be  a  flexible,  elastic  slab,  and 
the  subgrade  to  act  as  a  series  of  independent  springs.  Hogg  (k)   extended 
this  solution  to  the  case  where  the  top  layer  is  again  considered  to  be  a 
flexible,  elastic  slab,  but  where  the  subgrade  is  treated  as  a  semi- 
infinite  elastic  foundation.  A  further  refinement  was  introduced  by 
Burmister  (5)  who  solved  the  problem  of  the  two-layer  system  by  applying 
the  governing  equations  of  elasticity  to  both  layers.  He  assumed  that 
each  layer  acts  as  a  continuous,  isotropic,  homogeneous,  linearly  elastic 
medium  infinite  in  horizontal  extent.  His  methods  were  later  extended  to 
cover  more  than  two  layers  and  a  wider  range  of  relative  thicknesses  and 
material  constants  (6,  7,  8,  9>  10). 

The  important  consideration  in  the  above  discussion  is  that  all  of 
these  methods  are  based  on  the  condition  of  ideally  elastic  materials, 
with  fixed  elastic  constants.  That  no  material  is  perfectly  elastic  is 
generally  agreed,  and,  of  course,  highway  materials  are  no  exception. 
Under  a  fast  moving  wheel  load,  at  freezing  temperatures,  nearly  complete 


elastic  behavior  may  be  observed.  However,  the  ultimate  behavior  of  a 
pavement,  and  particularly  of  an  asphaltic  concrete  pavement,  is  obviously 
nonelastic,  since,  over  a  number  of  years,  and  under  the  action  of  a  great 
many  wheel  loads,  the  individually  insignificant,  nonrecoverable  strains 
accumulate  and  are  felt  in  many  forms,  such  as  rutting  in  the  wheel  path, 
flow  in  the  direction  parallel  to  that  of  the  traffic,  or  simply  as  a 
brittle  crack. 

Use  of  Viscoelastic  Theory 

In  order  to  take  into  consideration  the  time  dependent  response  to 
loads  of  the  pavement  materials,  they  are  generally  assumed  to  be  linear- 
ly viscoelastic.  This  simplifies  their  representation  and  makes  analyses 
less  complex.  Stress  and  strain  for  linear  viscoelastic  materials  can  be 
related  by  either  differential  or  integral  linear  operators  (ll).  The 
differential  operator  form  of  the  stress -strain  law  is  most  commonly  used 
and  as  Biot  (12)  has  pointed  out  may  be  visualized  as  a  combination  of 
springs  and  dashpots . 

For  a  number  of  years  engineers  have  recognized  that  the  component 
materials  of  a  pavement  system  have  time  dependent  stress-strain  proper- 
ties due  to  consolidation  and  creep  (13,  Ik).     Layered  systems  were  not, 
however,  analyzed  using  time  dependent  (viscoelastic)  material  properties 
until  after  the  development  of  the  correspondence  principle  for  isotropic 
media  in  1955  by  Lee  (ll)  and  shortly  thereafter  the  extension  by  Biot 
(15)  to  include  anisotropic  media.  More  general  methods  of  viscoelastic 
stress  analysis  have  been  described  by  Lee  (16)  and  Lee  and  Rogers  (17). 

Viscoelastic  layered  systems  have  been  solved  for  only  two  layers 
characterized  by  idealized  stress-strain  material  properties  such  as  those 


represented  by  Maxwell,  Voigt,  and  four  element  viscoelastic  models. 
Most  of  these  solutions  consider  a  Winkler-type  foundation  (l8,  19,  20, 
21),  or  a  system  of  linear  elements  with  only  limited  interaction  (22,  23, 

Monismith  and  Secor  (25)  attempted  to  validate  viscoelastic  theory 
for  a  simple  layered  system  by  comparing  theoretical  predictions  of  sur- 
face deflection  with  experimental  values.  The  system  consisted  of  a 
square  slab  of  asphaltic  concrete  which  rested  on  a  bed  of  closely  spaced 
springs.  It  was  observed  that  the  agreement  between  the  experimental  and 
theoretical  deflections  was  not  very  good,  especially  at  higher  test 
temperatures . 

Perloff  and  Moavenzadeh  (26)  presented  an  analysis  for  deflection 
due  to  a  point  load  moving  across  the  surface  of  a  semi-infinite  homoge- 
neous linear  viscoelastic  medium  assuming  the  response  of  the  material  to 
deviatoric  stress  to  be  that  of  a  Kelvin  model.  Using  an  N-element  Kelvin 
model,  Barksdale  and  Leonards  (27)  developed  a  method  for  viscoelastic 
analysis  to  investigate  pavement  systems  subjected  to  repeated  loadings 
and  found  that  this  method  gave  a  realistic  first  approximation  to  the 
actual  behavior  of  a  bituminous  pavement  system.  Harr  (28)  used  a  two 
element  model  to  show  the  influence  of  vehicle  speed  on  pavement  deflec- 
tion. 

Other  models  ranging  from  very  simple  ones  to  infinite  series  of 
springs  and  dashpots  have  been  utilized  in  attempts  to  explain  the  time 
dependent  deformations  in  pavements.  To  date,  however,  no  single  visco- 
elastic model  that  would  permit  accurate  prediction  of  pavement  perfor- 
mance has  emerged. 


Material  Coefficients  for  Transverse  Anisotropy 
The  design  of  a  pavement  system,  like  that  of  other  engineering 
structures,  must  start  with  an  understanding  of  the  materials  involved. 
In  both  the  elastic  and  viscoelastic  approaches  the  asphaltic  concrete 
has  been  assumed  to  be  homogeneous  and  isotropic.  However,  evidence 
gathered  by  researchers  indicates  that  bituminous  mixtures  as  used  in 
pavements  do  not  conform  to  the  behavior  predicted  by  the  assumption  of 
isotropy.  Roscoe  and  Schofield  (29)  and  others  (30,  31)  have  promoted 
concepts  of  anisotropic  behavior  in  soils  and  bituminous  mixtures.  Nevitt 
(32)  and  Gaudette  (33)  have  shown  that  particle  orientation  takes  place 
during  compaction  of  bituminous  mixtures.  In  the  pavement,  aggregate  par- 
ticles tend  to  arrange  themselves  in  positions  with  the  long  axis  hori- 
zontal. This  preferred  orientation  is  a  cause  of  transversely  anisotropic 
behavior . 

Hearmon  (3*0  presented  a  symmetric  matrix  established  with  the  aid 
of  strain-energy  functions  for  a  homogeneous,  elastic,  transversely  aniso- 
tropic material  which  conserves  strain  energy.  This  matrix  has  five  in- 
dependent coefficients  which  describe  the  stress-deformation  behavior  of 
such  a  material. 

Recognizing  the  fact  that  strain  energy  is  not  conserved  in  bitumi- 
nous mixtures  whose  stress-deformation  behavior  is  time  dependent, 
Busching  (35)  derived  a  matrix  of  material  coefficients  without  the  use 
of  strain-energy  functions.  The  resulting  matrix  is  nonsymmetric  and  has 
six  independent  coefficients  that  describe  the  behavior  of  a  homogeneous, 
transversely  anisotropic,  nonconservative  material  subject  to  small 
strains.  These  coefficients  are  functions  of  time  and  temperature.  The 


stress-strain  constitutive  relations  in  terms  of  these  six  material  coef- 
ficients are  as  follows : 
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€2    =    -IT       ST  a2    +       ET  (2) 


c      ,    Vl  .     ^3l°2    +      1  _ 

€3    '    ~  +    ~    +    E^°3  (3) 


%      -     G^%  00 


€5     =     G^a5  (5) 


€6    =    E^  (1  -  Ma.)  *6  <6) 


The  six  material  coefficients  are  easily  determined  for  a  given  ma- 
terial at  a  given  temperature  using  uniaxial  compression  and  shear  tests 
as  a  function  of  time.  Hence,  Eqs.  1  through  6  can  be  used  to  predict 
the  displacements  at  any  point  knowing  the  stresses  at  that  point.  How- 
ever, from  a  practical  standpoint  there  are  limitations  to  the  use  of 
these  equations.  Busching  evaluated  the  coefficients  over  a  period  of 
only  30  seconds  because  it  was  not  possible  to  evaluate  the  shear  modulus 
G.  over  longer  duration,  though  it  was  possible  with  the  other  coeffi- 
cients. Besides,  these  equations  cannot  be  applied  for  dynamic  loadings. 
Hence,  in  a  time -dependent  analysis,  such  as  in  the  case  of  bituminous 


mixtures,  Eqs.  1  through  6  have  limited  application. 

Transfer  Function  as  a  Material  Characteristic 
It  is  obvious  from  the  above  discussions  that  none  of  the  current 
theoretical  methods  adequately  describe  the  time-dependent  response  of 
bituminous  mixtures  and  that  there  exists  a  need  to  find  a  parameter  or 
parameters  or  a  function  which  can  satisfactorily  serve  the  purpose.  Such 
a  function  should  necessarily  be  time  dependent  and  should  be  able  to  be 
used  as  the  material  characteristic  for  both  static  and  dynamic  loads. 

A  study  of  the  literature  reveals  that  transfer  functions  are  fre- 
quently used  to  characterize  a  dynamic  system  in  electrical  and  mechanical 
engineering  problems  (36,  37,  38>  39)-  The  transfer  function  serves  as 
the  link  between  the  time -dependent  input  and  output  of  the  system.  A 
simple  example  is  that  of  a  mechanical  system  consisting  of  a  spring  with 
a  spring  constant  K,  a  dashpot  with  a  viscosity  coefficient  B,  and  a  mass 
M,  subjected  to  a  forcing  function  f(t)  which  is  the  input  to  the  system 
(39).  The  differential  equation  of  motion  for  this  system  is  given  by 

Mx"(t)  +  Bx(t)  +  Kx(t)  =  f(t)  (7) 

where  x(t)  -  output  of  the  system  which  is  the  displacement  from  the 
position  of  stable  equilibrium; 
x(t)  =  velocity  of  the  system  with  respect  to  the  inertial  refer- 
ence. 
x(t)  =  acceleration  of  the  system. 
The  transfer  function  for  this  system  is  obtained  from  the  ratio  of  the 
Laplace  transform  of  the  output  x(s)  to  that  of  the  input  f(s),  as 

G(s)=JJkL  (8) 

f(s)  . 
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Laplace  transform  of  Eq.  7,  neglecting  initial  conditions,  yields  the  re- 
lationship between  x(s)  and  f(s)  as 

(Ms2  +  Bs  +  K)  x(s)  =  f(s)  (9) 

where  s  is  a  complex  variable.  From  Eq.  Q,  the  transfer  function  for  the 
system  is  seen  to  be 

8(a)  =il£L=  _1 (9.1) 

f(s)     Ms  +  Bs  +  K 

Simple  as  well  as  complicated  electrical  networks  are  also  analyzed 
through  the  use  of  transfer  functions  to  study  their  frequency  responses 
(36,  39)'  Also,  transfer  functions  have  been  successfully  used  to  analyze 
such  diverse  dynamical  problems  as  the  rolling  motion  of  an  airplane  (kO) 
on  one  hand,  and  the  shock  and  vibration  in  mechanical  systems  (Ul)  on  the 
other . 

In  the  light  of  the  above  discussion,  the  transfer  function  appears 
to  have  the  potentialities  of  the  function  that  is  needed  to  study  the 
dynamic  response  of  viscoelastic  materials.  The  present  investigation 
describes  a  method  of  obtaining  the  transfer  function  for  a  given  bitumi- 
nous concrete  from  laboratory  tests  and  studies  the  extent  to  which  this 
function  can  represent  the  viscoelastic  characteristics  of  the  material. 
Effects  of  mix  type,  temperature  and  anisotropy  on  the  transfer  function 
were  studied  with  a  view  to  providing  insight  into  this  function  which, 
it  is  hoped,  may  facilitate  the  development  of  more  refined  methods  for 
the  structural  design  of  pavements. 


TRANSFER  FUNCTIONS  FOR  VTSCOELA.STIC  MATERIALS 

The  possibility  of  applying  the  concept  of  transfer  functions  for 
viscoelastic  materials  was  put  forth  in  the  previous  section.  The  means 
of  achieving  this  and  the  assumptions  involved  are  discussed  in  this 
section. 

Concepts  of  a  Mechanical  System 

The  applied  stresses  due  to  wheel  loads,  dynamic  or  static,  and  the 
resulting  displacements  are  the  subjects  of  interest  in  the  case  of  bi- 
tuminous mixtures  placed  in  a  pavement  structure.  This  is  analogous  to  a 
mechanical  system  subjected  to  a  force  input  and  a  displacement  output. 
Inasmuch  as  the  dynamical  behavior  of  such  a  mechanical  system  can  be 
easily  analyzed  through  transfer  functions,  it  is  logical  to  consider  the 
bituminous  mixtures  as  mechanical  systems  for  the  purpose  of  stress- 
d  i  splac  ement  analy s  i  s . 

The  three  fundamental  dynamical  properties  of  any  mechanical  system 
are  inertia,  damping  and  restoring  force.  These  three  are  represented  by 
the  three  basic  parameters  of  systems,  mass,  damping  function,  and  spring 
function  respectively.   If  the  parameters  are  isolated,  then  the  system 
is  called  a  lumped-parameter  system.  To  say  that  the  parameters  are  iso- 
lated is  to  mean  that  all  of  the  mass  of  the  system  is  lumped  into  rigid 
bodies;  all  of  the  damping  property  of  the  system  is  lumped  into  dampers; 
and  the  entire  spring-like  property  of  the  system  is  lumped  into  springs. 
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However,  in  actuality,  a  physically-real  spring  has  mass  and  is  also  a 
damper;  a  physical  damper  has  mass  and  some  spring-like  properties;  and  a 
physically  real  body  is  deformable  and  therefore  has  spring  and  damping 
properties.  To  say  that  a  system  is  a  lumped  parameter  one  is  therefore 
to  state  an  approximation.  The  dynamical  behavior  of  such  a  system  is 
described  by  a  second  order  linear  differential  equation. 

On  the  other  hand,  a  mechanical  system  in  which  the  basic  parameters 
are  not  isolated  is  called  a  continuous  or  distributed-parameter  system. 
The  dynamical  response  of  such  a  system  is  analyzed  through  a  differential 
equation  of  order  higher  than  two  which  represent  the  system  behavior. 

Asphaltic  Concrete  as  a  Continuous  System 
It  has  been  seen  from  the  discussion  under  "Use  of  Viscoelastic 
Theory"  that  attempts  to  define  the  response  of  bituminous  mixtures  by  the 
use  of  spring-dashpot  models,  that  is,  by  assuming  them  to  be  lumped  param- 
eter systems,  have  not  yielded  satisfactory  results.  Hence,  in  this  in- 
vestigation the  asphaltic  concrete  has  been  analyzed  as  a  continuous  me- 
chanical system. 

Depending  on  the  degree  of  damping,  such  a  system  can  be  underdamped, 
critically  damped  or  overdamped.  The  nature  of  the  system  can  easily  be 
recognized  through  the  amplification-frequency  plot  in  a  frequency  spec- 
trum analysis  (^l).  A  distinct  peak  in  the  spectrum  indicates  an  under- 
damped  system  whereas  the  absence  of  a  peak  and  continuous  fall  of  the 
curve  shows  an  overdamped  system. 

The  frequency  spectrums  for  the  asphaltic  concrete  specimens  tested 
in  this  investigation  indicate  an  overdamped  system. 


11 

Analysis  of  Mechanical  Systems  Through  Transfer  Functions 
It  has  been  pointed  out  in  the  section  on  "Transfer  Function  as  a 
Material  Characteristic"  that  transfer  functions  have  been  used  to  ana- 
lyse the  dynamical  response  of  mechanical  systems.  In  those  cases  where 
the  system  behavior  is  known,  the  transfer  function  is  easily  obtained 
from  the  differential  equation  relating  the  output 'and  the  input  vari- 
ables. However,  for  systems  whose  behavior  is  unknown,  the  transfer  func- 
tion may  be  approximated  from  certain  types  of  input-output  data  obtained 
from  laboratory  tests  (36,  39,  kl) .     This  method  of  obtaining  the  transfer 
function  is  followed  in  the  case  of  asphaltic  concrete  specimens  in  this 
investigation  in  the  absence  of  any  prior  knowledge  of  their  dynamical 
behavior . 

Further,  it  is  assumed  that  asphaltic  concrete  behaves  as  a  linear 
system.  Although  it  has  been  recognized  that  the  behavior  of  this  materi- 
al is  nonlinear  particularly  at  higher  stress  levels  and  at  higher  tem- 
peratures, the  investigations  of  Busching,  Goetz  and  Harr  (^2)  and  others 
{k3,   kk)   have  established  the  validity  of  this  assumption  for  small  dis- 
placements normally  encountered  in  bituminous  mixtures . 
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RELATED  THEORY 

It  has  been  pointed  out  in  the  previous  section  that  asphaltic  con- 
crete  can  be  considered  as  a  continuous  mechanical  system  in  order  to 
analyze  its  time -dependent  response  to  an  applied  load  and  that  transfer 
functions  can  be  used  to  characterize  this  response.  The  basic  theory  be- 
hind the  transfer  functions  is  presented  in  this  section. 

Concept  of  Transfer  Function 
The  ratio  of  an  operational  output  (Laplace  transform  of  the  output) 
of  a  dynamical  system  to  the  operational  input  (Laplace  transform  of  the 
input)  is  called  the  transfer  function  between  the  operational  input  and 
its  corresponding  output.  The  transfer  function  is  shown  schematically 
in  Fig.  1.   It  is  by  definition  a  function  of  the  complex  variable  s. 


Operational 


-*■ 


input 


Transfer 
function 


Operational 

*- 


output 


Fig.  1.   Block  diagram  representation 
of  the  transfer  function. 
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If  a  rigid  body  is  subjected  to  a  dynamical  force  f(t),  then  the 
force  and  the  resulting  displacement  x(t)  can  be  considered  as  the  input 
and  output  respectively  for  the  dynamical  system.  The  transfer  function 
between  the  operational  force  and  the  operational  displacement  is  given  by 

g(b)  =  3sl  (10) 

f(s) 

where       G(s)  =  the  transfer  function 

x(s)  =  Laplace  transform  of  the  output 

f(s)  =  Laplace  transform  of  the  input. 
Solving  for  x(s), 

x(s)  =  G(s)f(s)  (11) 

The  inverse  transform  of  x(s)  is  x(t)  and  is 

x(t)  =  £e'(G(s)f(s)}  (12) 

Eq.  12  shows  that  once  the  transfer  function  G(s)  is  known  for  any  system, 
the  displacement  x(t)  can  be  evaluated  for  that  system  for  another  given 
input  force  f (t) . 

In  general,  the  transfer  function  G(s)  is  a  ratio  of  two  polynomials 
in  s,  such  as 

r(   >    Ms)    (s  -  aj)(s  -  ap  ...  (s  -  a^) 


D(s7  =   (S  -  b£)(8  -  b^  ...  (S  -  b^) 


where  a.'  's  and  b.'  *s  are  the  roots  or  zeros  of  the  numerator  N(s)  and 
i        l 

denominator  D(s)  of  G(s)  respectively,  b  '  's  are  also  called  the  poles 


Ik 

of  G(s),  since  the  transfer  function  becomes  infinite  when  evaluated  at  a 
zero  of  its  denominator.  That  is, 


G(S)  .-  cc  (lit) 

s— b. 


Substituting  for  G(s)  in  Eq.  11,  x(s)  can  be  rewritten 


f(s)(s  -  ajMs  -  a')...(s  -  a') 
X(S)  =    («-b')(,  -4)...(s-b;)  («) 


Eq.  15  can  be  expanded  in  m  partial  fractions  for  the  m  roots  of  D(s)  as 

*v--zhrrj  *  Tr^Ep*  —  +Tr?VT    (16) 

l  c.  m 

which  is  an  identity  in  s  where  C',  C',  ,  C'  are  constants  independent 

of  s.  Using  the  Laplace  inverse  transform,  the  final  solution  for  x(t)  is 
obtained  from  Eq.  16  (i+l)  as 


b'  t      b'  t  b'  t 

x(t)  =  c'e  L     +  C'e         +  ...  +  C'e  m  (17) 

v  '  1       2  m  v  " 


From  the  above  discussion  it  is  obvious  that  once  the  transfer  func- 
tion for  a  system  is  known,  the  displacement  x(t)  of  the  system  for  a 
given  force  input  f(t)  can  be  evaluated  theoretically.  The  solutions  for 
displacements  for  two  specific  cases  of  interest,  namely,  a  sinusoidal 
load  input  and  a  step  function  input,  are  presented  below  for  a  system 
whose  transfer  function  is  given  by  G(s). 
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Solution  for  a  Sinusoidal  Input 
Let  the  force  input  for  the  system  whose  transfer  function  is  given 
by  G(s)  be 

f(t)  =  f  sin  tot,  (16) 

where  f  is  the  magnitude  of  the  sine  function  and  co  is  the  time  frequen- 
cy. The  operational  displacement  x(s)  of  the  system  is  given  by  Eq.  11, 

x(s)  =  G(s)f(s)  (11) 

where 

f(s)  =<£{f(t)} 

=  c£{f0  sin  "*) 


f  CO 

o 


2  .   2 
s  +  to 


(19) 


Substituting  Eq.  19  in  Eq.  11  and  solving  for  x(t)  by  using  the  method  of 
partial  fractions, 


x(t)  =0C1{x(s)} 


-IfG(s)  tn(D 


12.2 

^  s  +  to 

=  f  |G(Jcd)  J  sin  (tot  +  ♦),    (20) 

where  4>  is  the  phase  angle  between  the  input  and  the  output,  and  t  repre- 
sents time . 
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Solution  for  a  Step  Function  Input 
A  step  function  input  is  defined  as 

f(t)  =0      f or  t  <  0 


f(t)  =  f 


for  t  ^  0 


(21) 


where  f  is  a  constant  with  time  as  shown  in  Fig.  2.  In  view  of  this 
o  < 

definition,  for  large  times,   a  static  compressive  force  of  magnitude  f 
applied  to  a  system  can  be  considered  as  a  step  function  input. 


Fig.  2.  A  step  function  of  time. 

The  Laplace  transform  of  the  step  function  f (t)  is  given  by 

f(s)  =<£{f(t)} 


(22) 


As  cited  previously,  the  operational  displacement  of  a  system  subjected 
to  any  force  input  is  given  by  Eq.  11, 

x(s)  =  G(s)f(s)  (11) 
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Substituting  Eq.  22  into  Eq.  11, 


f 
x(s)   =  G(s)-£-  (23) 


s 


The  inverse  transform  of  x(s)  is  x(t)  which  is  obtained  from 

-iff  G(sh 
I(t)  .£      ^ 

Assuming  G(s)  of  form 


(21+) 


A(s  -  ax)(s  -  a2)(s  -  a  ) 
G(S)  =  (s  -  blKs  -  b2)(s  -  b3)(s  -  bu)     (25) 


as  is  the  case  in  the  present  investigation,  where  A  is  a  constant,  a.  's 
and  b.  's  are  zeros  of  the  numerator  and  denominator  respectively,  x(t)  is 
obtained  from  Eq.  2k,   as 

_l  ;A  fQ(s  -  ax)(s  -  a2)(s  -  a  )     ^ 
X(t)  =X  \(s  -  bl)(s  -  b2)(s  -  b3)(s"-  bu)s|    (2b) 

Eq.  26  can  be  solved  using  the  method  of  partial  fractions  which  yields 


b  t     b  t     b  t     b,t 
x(t)  -  C^e    +  C2e    +  c3e    +  cl+e    +  c5         (2?) 


where  C.  's  are  constants  independent  of  t. 

Determination  of  the  Transfer  Function 
The  experimental  basis  to  obtain  the  transfer  function  is  the  fre- 
quency spectrum.  A  sinusoidal  input  of  known  magnitude  and  frequency  is 
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applied  to  the  system  whose  transfer  function  is  sought.  The  magnitude 
of  the  steady-state  output  and  the  phase  angle  between  the  output  and  the 
input  are  then  measured.  It  is  necessary  to  make  certain  that  the  output 
is  only  the  steady-state  component  (after  the  transient  component  has 
been  assumed  to  have  damped  out)  by  measuring  its  frequency  and  ascertain- 
ing that  it  is  the  same  as  that  of  the  input.  In  a  frequency  spectrum 
analysis,  interest  is  in  the  ratio  of  magnitude  and  not  in  the  value  of 
each  of  the  magnitudes  of  the  input  and  steady-state  output. 


Input 


Output 


amplification  =  a2/ai 
phase  angle  =  ojt 


Fig.  3.   Experimental  magnitude  and  phase  measurements 


The  ratio  of  the  magnitude  of  the  output  to  that  of  the  input  is 
called  the  amplification  or  gain.  The  amplification  and  phase  angle  are 
therefore  experimentally  determined  at  a  given  frequency.  The  frequency 
of  the  input  is  then  changed  to  another  value  and  again  the  amplification 
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and  phase  angle  are  determined.  The  procedure  continues  from  a  very  low 
frequency  (approximately  zero)  to  that  frequency  at  which  the  output  mag- 
nitude is  approximately  zero.  The  plot  of  the  amplification  against  fre- 
quency is  called  the  frequency  spectrum,  and  for  convenience  the  magni- 
tude  of  the  amplification  is  plotted  in  decibels  .  Fig.  3  shows  the  in- 
put  and  output  traces  for  a  given  frequency. 

Once  the  amplification  is  plotted  against  log  of  frequency,  the  curve 
may  be  approximated  with  a  series  of  connected  straight  lines  as  shown  in 
Fig.  h.     These  straight  lines  are  selected  so  as  to  form  asymptotes  to  the 
curve.  Although  straight  lines  of  any  slope  nay  be  used,  the  slopes  of 


Amplification,  decibels 


^L 


«  0 


k.  =  slope,  decibels/octave 


cd.  =  corner  frequency,  rad/sec 


Fig.  k.     Asymptotic  approximation  of  a  frequency  spectrum. 


Rdb  "  20  lo^io  R 
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the  asymptotes  are  generally  taken  to  be  six  decibels  per  octave  for 
greater  ease  of  determination.  The  closer  the  approximation  to  the  actual 
frequency  spectrum  curve,  the  better  will  be  the  approximation  of  the 
transfer  function  being  sought  (Ul). 

The  intersections  of  the  selected  approximating  asymptotes  determine 
the  corner  frequencies.   In  Fig.  k   these  are  six:  ax.    =0,  cop,  o5_,  m  ,  ov 
and  o^.  The  slopes  of  the  asymptotes  are  given  by  k_  ,  k  ,  k  ,  k,  ,  k  and 
k/-  decibels  per  octave.   Once  the  corner  frequencies  and  the  slopes  of  the 
approximating  asymptotes  are  obtained,  the  modified  transfer  function 
which  describes  the  frequency  spectrum  is  written  (Ul)  as 

kl  k2'h  Vk2  Vkn-l 

G(jcd)   =  A(jco  +  a^)        (ja>  +  a^)  (ja>  +  ol)  .  . .    ( ja>  +  a>n)  (28) 

The  modified  transfer  function  is  the  vector  amplification  or  simply  the 
ratio  of  the  magnitude  of  the  output  to  that  of  the  input  and  is  a  func- 
tion of  the  input  frequency.  By  substituting  s  =  jco  in  Eq.  28,  the 
transfer  function  for  the  system  is  obtained 


*1       Vkl       Yk2  Vkn_l 

G(s)  *A(s  +<dl)~5(s  +ov,)     (s+o^)     ...(s+o>n)        (29) 


The  constant  A  can  be  determined  from  any  one  experimental  frequency  in 
the  frequency  spectrum  by  calculating  the  absolute  value  of  the  transfer 
function  from  Eq.  28  at  that  frequency  and  equating  it  to  the  correspond- 
ing amplification  value  in  the  frequency  spectrum. 


4f 

An  octave  is  defined  as  the  frequency  range  for  which  the  ratio  of 

the  upper  bound  frequency  to  the  lower  bound  frequency  is  2  to  1; 

that  is,  a).  <  cu  <  2to.  is  an  octave. 
'  i        1 
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A  typical  example  of  determining  the  transfer  function  and  the  con- 
stant A  from  a  frequency  spectrum  obtained  for  a  bituminous  concrete  speci- 
men is  given  in  the  Appendix. 

It  must  be  noted  that  this  method  of  obtaining  the  approximate  trans- 
fer function  from  the  frequency  spectrum  is  applicable  only  for  an  over- 
damped  or  critical ly  damped  system.   It  may  be  used  in  the  present  inves- 
tigation because  bituminous  mixes  have  been  observed  to  behave  as  over- 
damped  systems. 
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EXPERIMENTAL  INVESTIGATION 

The  experimental  phase  of  this  study  had  as  an  objective  the  develop- 
ment and  use  of  accurate  techniques  for  obtaining  all  of  the  data  required 
in  the  determination  of  the  transfer  function  from  a  frequency  spectrum. 
The  core  of  the  technique  lies  in  devising  a  system  which  can  apply  a  si- 
nusoidal force  input  of  a  desired  magnitude  and  frequency  and  which  can 
simultaneously  measure  the  displacement  output. 

An  MTS  electronic  function  generator  coupled  with  a  loading  frame 
which  was  fitted  with  an  electronically  controlled  hydraulic  actuator  was 
used  for  this  purpose.   See  Figs.  5  and  6.  The  loading  frame  was  in- 
stalled in  a  constant  temperature  room  which  facilitated  carrying  out  the 
tests  at  the  desired  temperatures .  A  two-channel  Brush  recorder  which 
was  attached  to  the  function  generator  recorded  the  input  and  output  of 
the  test  specimen  simultaneously. 

Dynamic  Tests 

The  function  generator  was  used  to  apply  sinusoidal  force  inputs  to 
the  asphaltic  concrete  specimen  by  suitable  manipulation  of  the  controls . 
The  magnitudes  of  the  sine  force  varied  from  2.5  to  15  lbs.  and  the  fre- 
quencies varied  from  0.0016  to  1.6  cycles  per  second. 

Since  in  a  sine  test  the  specimen  is  alternatively  in  compression 
and  in  tension,  the  method  of  mounting  the  specimen  in  the  test  assembly 
was  of  paramount  importance.  Two  aluminum  plates,  2  in.  x  2  in.  x  3/8  in. 
in  size,  were  glued  to  the  specimen  top  and  bottom  respectively,  using 
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commercially  available  Eastman  910  Adhesive.  See  Fig.  7.  The  plates,  in 
turn,  were  connected  to  the  top  and  bottom  loading  platens  by  screws  and 
nuts.  The  top  platen  was  rigidly  fixed  to  the  moving  ram  of  the  hydraulic 
actuator  while  the  bottom  platen  was  fixed  to  the  loading  deck.  See  Fig. 
8.  Remotely  controlled  and  commanded  by  the  electronic  console  of  the 
function  generator,  the  hydraulic  ram  moves  up  and  down  at  the  pre-set 
sine  force  and  frequency.  Thus  the  specimen  was  subjected  to  the  desired 
sinusoidal  force  input  which  was  recorded  by  the  left  channel  of  the  Brush 
recorder.  The  corresponding  displacement  of  the  specimen  was  sensed  by 
the  LVDT  in  the  actuator  and  was  simultaneously  recorded  as  the  output  in 
the  right  channel  of  the  recorder.  Fig.  9  shows  the  two  traces  thus  ob- 
tained in  a  typical  dynamic  test. 

Static  Tests 
The  same  test  assembly  which  was  used  for  the  dynamic  tests  was  used 
for  the  static  tests.  When  the  function  generator  is  set  for  a  combina- 
tion of  a  ramp  function  and  a  high  frequency,  the  actuator  applies  a  stat- 
ic compressive  force  of  desired  magnitude.  As  before,  the  left  and  right 
channels  recorded  the  input  and  output  respectively.  The  two  traces  of  a 
typical  static  compression  test  are  illustrated  in  Fig.  10. 

Materials  and  Preparation  of  Specimens 
Two  different  gradings,  one  with  a  maximum  size  of  the  No.  k   sieve 
and  the  other  the  3/8  inch  sieve  were  used  in  this  investigation.  The 
bituminous  mixes  prepared  on  the  basis  of  these  gradings  were  designated 
as  Mix-1  and  Mix-2,  respectively.  Table  1  shows  the  sieve  analysis  and 
origin  of  the  aggregates  used.  The  grading  curves  are  given  in  Fig.  11. 
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Fig.  7.  A  Typical  Test  Specimen 


Fig.  8.   Specimen  in  Position  for  Testing 
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Fig-  9-  Typical    traces     of    a    sinusoidal    loading    test 
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Fig.  10.  Typical    traces     of   a    static    compression    test 


29 


TABLE 
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Sieve 

Analysis  of 
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No. 
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No. 

6 

No. 
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A  60-70  penetration  asphalt  cement  having  the  properties  shown  in  Table  2 
was  used.  After  being  blended  in  2200  gram  batches,  the  aggregate  was 
placed  in  an  oven  and  heated  to  325°F.  The  asphalt  was  heated  to  325°F 
and  then  mixed  with  the  heated  aggregate  in  a  mechanical  mixer  for  two 
minutes.  The  weight  of  asphalt  cement  mixed  with  the  aggregate  was  5.0 
and  3.5  percent  by  weight  of  the  aggregate  for  Mix-1  and  Mix-2,  respec- 
tively. After  mixing  was  completed,  each  batch  was  put  in  flat  pans  and 
placed  in  a  forced-draft  oven  to  cure  for  2k   hours  at  lU0°F. 

Compaction  of  the  bituminous  mixture  was  started  after  the  mixture 
was  reheated  to  325°F.  The  heated  mixture  was  placed  in  two  layers  in 
the  heated  beam  mold  of  internal  dimensions  2  in.  x  2  l/2  in.  x  12  in. 
Each  layer  was  rodded  25  times  with  a  l/2-inch  diameter,  round-end  tamping 
rod.  The  mixture  in  the  beam  mold  was  compacted  with  the  kneading  com- 
pactor by  the  use  of  a  special  rectangular  tamping  foot  of  size  2  in.  x 
2  l/2  in.  The  mixture  was  compacted  with  the  tamping  foot  by  moving  the 
mold  back  and  forth  its  entire  length  under  the  tamping  foot  four  times, 
using  eight  tamps  during  each  one-way  traverse  for  each  of  four  different 
compaction  pressures.  Foot  pressures  used  in  compaction  were  20,  35»  100 
and  120  psi  in  that  order.  After  the  beam  was  compacted  on  its  top  side 
according  to  this  procedure,  the  mold  and  beam  were  inverted  and  the  pro- 
cedure was  repeated.  Hence,  both  the  top  and  the  bottom  of  the  beam  re- 
ceived the  same  compactive  effort. 

After  the  beam  was  compacted,  it  was  removed  from  the  mold  and  placed 
in  a  forced-draft  oven  and  cured  for  2k   hours  at  a  temperature  of  lU0°F. 
After  the  curing  period,  the  beam  was  cooled  to  room  temperature. 

The  compacted  beam  was  cut  into  specimens  by  use  of  a  masonry  saw. 
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TABLE  2 
Results  of  Tests  on  Asphalt  Cement 


Specific  Gravity  at  77°F 1.036 

Softening  Point,  Ring  and"  Ball,  °F 12U 

Ductility  at  77°F,  5  cm/min,  cm 100+ 

Penetration,  100  grams,  5  sec,  77°F 63 

Penetration,  200  grams,  60  sec,  32CF 17 

Loss  on  Heating,  50  grams,  5hr.,  325°F,  # 0.01 

Penetration  of  Residue,  %   of  Original  89 

Flash  Point,  Cleveland  Open  Cup,  °F 595 

Solubility  in  CCl^,  % 99-84 
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One-inch  pieces  were  cut  from  both  ends  of  the  beams  and  wasted.  This 
was  done  because  it  was  recognized  that  confinement  at  the  ends  of  the 
beam  mold  might  cause  compaction  conditions  to  be  atypical  at  these  lo- 
cations. Specimens  were  cut  along  the  x,  y,  and  z  axes  of  the  beam,  as 
shown  in  Fig.  12,  approximately  0.05  inch  oversize.  Final  shaping  and 
polishing  were  accomplished  on  a  lapping  wheel  using  100  and  800  mesh 
grinding  compounds  successively.  The  dimensions  of  the  test  specimens 
were  controlled  to  0.01  inch  by  utilizing  the  jig  shown  in  Fig.  13,  and 
the  final  dimensions  were  1  in.  x  1  in.  x  2  in. 


-,  y 

f— *■   X 

1 

z 


Fig.  12.  Designation  of  axes  for  the  beam. 


Scope 
The  investigation  was  confined  to  the  two  bituminous  mixes  described 
in  the  previous  section  and  to  three  temperatures,  namely,  70,  80  and  90°F. 
For  each  mix  three  different  specimens  cut  along  the  x,  y  and  z  directions 
of  the  beam  were  used.  The  specimens  were  designated  as  x-1,  y-1,  z-1, 
x-2,  y-2  and  z-2,  the  letter  denoting  the  axis  along  which  it  was  cut  and 


Fig.  13.   Jig  for  Specimen  Alignment 
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the  numeral  the  mix. 

In  order  to  investigate  the  effect  of  size  of  specimen  on  the  trans- 
fer function,  three  different-sized  specimens  from  Mix-1  were  tested  at 
one  temperature,  namely,  90°F-  Keeping  the  height-width  ratio  constant, 
the  following  three  sizes  were  used:  3A  i-n-  x  3/^  in.  x  1  l/2  in.,  1  in. 
x  1  in.  x  2  in.  and  1  l/k   in.  x  1  l/k   in.  x  2  1/2  in. 
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RESULTS  OF  EXPERIMENTS 

The  results  of  all  the  dynamic  and  static  tests  at  70,  80  and  90°F 
are  summarized  in  Tables  3  through  10.  The  transfer  functions  derived  for 
each  specimen  at  these  temperatures  are  shown  in  Tables  3,  k   and  5.  The 
calculated  displacements  for  each  specimen  for  three  different  sinusoidal 
load  inputs  and  constant  load  inputs,  respectively,  and  the  corresponding 
measured  displacements  are  shown  in  Tables  6,  7  and  8. 

The  phase  angles  obtained  in  the  dynamic  tests  for  each  specimen  are 
shown  separately  in  Table  9  with  respect  to  the  frequency  of  the  test. 
The  same  data  are  shown  graphically  in  Figs .  Ik   through  19 • 

The  absolute  value  of  G(s)  versus  frequency  curves  for  the  six  speci- 
mens from  the  two  different  mixes  tested  at  the  three  different  tempera- 
tures are  shown  in  Fig.  22. 
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TABLE  3 
Transfer  Functions  at  70°F 


Specimen  Transfer  Function 

r(^   -     (0»03?)(s  +  0.03l)(s  +  Q.l8)(s  +  1.7) 
^S)   ~  (s  +  0.02)(s  +  0.1)(s  +  l)(s  +  10) 


G(S)  =  (0;02^)(s  +0.02)(s  |0.l8)(s  +  1.6) 
^S)  (s  +  0.0lU)(s  +  0.1)(s  +  l)(s  +  10) 

r(*)         (0.03^8)(s  +  0.03)(s  +  0.1Q(s  +  2.7) 

V  ;  '   (s  +  0.02)(s  +  0.25)(s  +  2)(s  +  10) 

_/  v    (0.032)(s  +  0.025) (s  +  0.l6)(s  +  1.8) 

G(S}  =   (s  +  0.02)(s  +  0.lj(s  +  l)(s  +  1C) 

(    v    (0.03^6)(s  +  0.023)(s  +  Q.l6)(s  +   1.6) 

^s;  ~   (s  +  0.02)(s  +  0.1)(s  +  l)(s  +  10) 

,  ,    (O.Q3^6)(s  +  0.023)(s  +  Q.l6)(s  +  1.6) 

{jKS)  =   (s  +  0.02)(s  +  0.1)(s  +  l)(s  +  10) 
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z-l 


x-1 


y-i 


z-2 
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TABLE  h 
Transfer  Functions  at  80°F 


Specimen  Transfer  Function 

r(,\  -      (0-05)(s  +  0.023)(s  +  0.2)(s  +   2) 

^S)  ~      (s  +  0.02)(s  +  0.1)(s  +  l)(s  +  10) 

(    v    (0.0m15)(s  +  0.028)(s  +  0.17)(s  +  1.8) 

U^S;  "   (s  +  O.02)(s  +  0.1)(s  +  l)(s  +  10) 

nt    \  (0.0*0  (s  +  0.02^)(s  +  Q.l6)(s  +  2.3) 

G(S}  =  (s   +  0.02)(s  +  0.l)(s  +  l)(s  +  10) 

-,  v    (O.OUHs  +  0.025) (s  +  0.16) (s  +2.3) 

G^  =   (s  +  0.02)(s  +  0.l)(s  +  l)(s  +  10) 

„,  v    (0.0Ul5)(s  +  0.028)(s  +  0.17)(s  +  1.8) 

G(SJ  =   (s  +  0.02) (b  +  0.l)(s  +  l)(s  +  10) 

G(S)  _  (0-0^15)(s  +  0.028)(s  +  0.17)(s  +  1.8) 

{jKS}  "  (s  +  0.02)(s  +  0.l)(s  +  l)(s  +  10) 


y-2 
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TABLE  5 
Transfer  Functions  at  90 °F 


Specimen  Transfer  Function 

,  v    (Q.l)(s  +  0.02)(s  +  0.2)(s  +  1.7) 

^s'  =  (s  +  0.015)(s  +  0.1)(s  +  l)(s  +  10) 

(    v    (0.082)(s  +  0.03)(s  4  Q.2)(s  +  1.6) 

U[S)  ~   (s  +  0.02)(s  +  0.1)(s  +  l)(s  +  10) 

rl    \  (Q.l)(s  +  0.02)(s  +  0.17)(s  +  2.0) 

UU;  =  (s+  0.015)(s  +  0.1)(s  +  l)(s  +  10) 

r(0\  _  (0.07)(s  +  0.03)(s  +  0.17)(s  +  2.5) 

^{S)  ~      (s  +  0.02)(s  +  0.1)(s  +  l)(s  +  10) 

/  v    (0.08$)(s  +  0.02)(s  +  0.2)(s  +  1.8) 

^S;  "   (s  +  0.015)(s  +  0.1)(s  +  l)(s  +  10) 

/  v    (0.085)(s  +  0.02)(s  +  0.2)(s  +  1.8) 

L,{S)  ~  (s  +  0.015)(s  +  0.1)(s  +  l)(s  +  10) 
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TABLE  9 
Phase  Angles 


Specimen 

♦, 

deg 

• 

Temp 

a1 

;  a>, 

rad/sec 

°F 

0.025 

0.05 

0.1 

0.25 

0, 

•  5 

l.C 

) 

2.5 

5-C 

) 

10. 

,0 

70 

z-1 

18.6 

21.5 

20. 

,6 

17. 

,2 

lit. 

•  3 

11. 

,14 

ll.lt 

10. 

,0 

5, 

.7 

x-1 

llt.it 

lU.lt 

11, 

•  5 

10. 

M 

8.6 

.  6. 

,9 

5.7 

5- 

.7 

y-i 

17.8 

23.0 

23. 

,0 

21, 

■  5 

in, 

,u 

11. 

.5 

10.0 

7. 

,2 

5, 

■  1 

z-2 

18.3 

21.5 

23. 

,0 

18.6 

lit, 

.It 

11, 

■  5 

11.5 

7. 

■  3 

5. 

.7 

x-2 

18.2 

21.5 

18, 

,2 

14. 

,4 

11, 

•  5 

ll. 

5 

8.6 

7, 

,2 

5. 

•  7 

y-2 

18.2 

21.5 

18, 

.2 

18, 

,2 

11, 

.5 

9. 

,2 

7.2 

7, 

,2 

5. 

.7 

80 

z-l 

21.8 

23-9 

23. 

,8 

21, 

■  5 

19 

.5 

17. 

.8 

15-5 

13. 

,2 

12 

.0 

x-1 

18.0 

20.2 

22, 

•  5 

21, 

.3 

18, 

.6 

16, 

.0 

13.0 

11, 

.2 

10 

.0 

y-i 

32.2 

33-U 

29. 

,8 

28, 

.7 

28, 

.7 

2l4, 

.1 

21.5 

20, 

.1 

17 

■J: 

z-2 

28.7 

Uo.o 

3U. 

•  5 

3U, 

.5 

28 

.7 

25, 

.8 

21.2 

lit. 

.14 

11 

M 

x-2 

28.7 

28.7 

3U 

.5 

28 

.7 

22 

.2 

25 

.2 

lU.lt 

lit, 

.14 

11 

.It 

y-2 

28.6 

31.5 

31 

.5 

25 

.1 

23 

.0 

23. 

.0 

21.5 

21 

.5 

17 

.2 

90 

z-l 

18.2 

23.0 

23 

.0 

18 

.2 

11 

•  5 

11 

.5 

8.6 

7 

.2 

5 

.7 

x-1 

18.2 

25.8 

23 

.0 

17 

.2 

111 

.14 

11 

■  5 

8.6 

7 

.2 

5 

.7 

y-i 

18.2 

2*4.5 

23 

.0 

17 

.2 

lit 

.U 

11 

.5 

8.6 

7 

.2 

5 

.7 

z-2 

21.5 

23-0 

23 

.0 

17 

.2 

lit 

.It 

11 

•  5 

9-7 

7 

.2 

5 

•  7 

x-2 

21.5 

28.6 

23 

.0 

17 

.2 

in 

.3 

11 

•  5 

8.6 

7 

.1 

5 

.7 

y-2 

21.5 

25.8 

23 

.0 

15 

•  5 

17 

.2 

11 

.5 

11.5 

7 

.2 

5 

.7 
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TABLE  10 
Results  of  a  Typical  Sinusoidal  Test 
Specimen:         z-1 
Test  Temperature :  80°F 


CO 

r ad/ sec 

al 

mv 

a2 
mv 

a 
R  =  — 
al 

Rdb 

t 
sec 

0 

rad      deg . 

0.025 
0.050 

1000 

II 

22.0 
19.5 

0.0220 
0.0195 

-33.1520 

-34.4240 

15.30 
8.32 

0.382 
0.416 

21.8 

23.9 

0.100 
0.250 
0.500 

II 
II 
II 

17.0 

11.0 

9-0 

0.0170 
0.0110 
0.0090 

-35.3920 
-39.1720 
-4 O.916O 

4.15 
1.50 
0.68 

0.415 

0.375 
0.340 

23.8 
21.5 
19.5 

1.000 
2.500 
5.000 

II 
II 
II 

8.0 
6.5 
5.0 

0.0080 
O.OO65 
0.0050 

-41.9380 
-43.7420 
-46.0200 

0.311 
0.110 

0.046 

0.311 
0.271 
0.230 

17.8 

15.5 
13.2 

10.000 

II 

k.o 

0.0040 

-47.958O 

0.021 

0.210 

12.0 

LEGEND 


O)    = 

input  frequency,  rad/sec 

al  - 

input  magnitude,  millivolts 

a?  = 

output  magnitude,  millivolts 

R 

a2/al 

Rdb  = 

R  expressed  in  decibels 

= 

20  log1QR 

time  lag  betv/een  input  and  output,  sec, 
phase  lag  between  input  and  output 
cot  radians . 
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DISCUSSION  OF  RESULTS 

The  doncept  of  transfer  function  for  a  given  system  and  the  methods 
of  predicting  the  displacements  of  the  system  for  a  sinusoidal  and  step 
function  input,  respectively,  discussed  in  the  previous  sections  were 
applied  to  the  analysis  of  the  test  data  obtained  for  the  two  bituminous 
concrete  mixtures.  The  discussion  of  the  results  are  presented  in  this 
section. 

Frequency  Response  of  Asphaltic  Concrete 
In  most  of  the  dynamic  tests  a  magnitude  of  10  lbs.  was  employed  for 
the  sinusoidal  load  input  for  frequencies  covering  three  decades,  varying 
from  0.01  radians  per  second  to  10.0  radians  per  second.  These  were  ob- 
served to  be,  respectively,  the  slowest  and  the  fastest  frequencies  to 
which  the  test  specimens  were  responsive.  The  input  magnitude  and  the 
output  displacement  were  recorded  in  pounds  and  inches,  respectively. 
Also,  both  of  them  were  recorded  in  millivolts. 

The  ratio  of  the  output  to  the  input  expressed  in  decibels  (also 
referred  to  as  amplification  or  gain)  plotted  against  log  frequency  gives 
the  desired  frequency  spectrum.  The  phase  lag  or  simply  the  phase  angle 
for  each  frequency  was  also  computed  from 

♦  =  cut  (30) 

where   *  =  the  phase  angle  in  radians 

<o  =  the  test  frequency  in  rad/sec,  and 

t  =  time  lag  between  input  and  output  in  seconds . 
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The  phase  angle  is  also  plotted  against  log  frequency.  The  two  plots  to- 
gether represent  the  frequency  response  of  the  test  specimen  completely. 

Table  10  shows  typical  test  results  for  a  dynamic  test  and  Fig.  20 
gives  the  frequency  response  for  the  same  test.  It  is  seen  that  the  dis- 
placement of  the  specimen  is  largest  at  the  slowest  frequency  and  is  least 
at  the  fastest  frequency.  Thus  the  gain  continuously  falls  with  increase 
in  frequency.  It  was  discussed  before  with  reference  to  damped  systems 
that  such  a  curve  indicates  the  overdamped  nature  of  the  test  specimen. 

The  phase  angle  increases  with  increase  in  frequency  in  the  first 
decade  of  test  and  then  decreases  so  that  a  bell-shaped  curve  results 
when  phase  angle  is  plotted  against  frequency.  The  peak  has,  in  general, 
been  observed  to  occur  in  the  region  of  0.05  to  0.1  rad/sec.   See  Figs. 
Ik   through  19.  These  results  substantiate  the  observations  of  Pagen  (U5) 
and  others  who  have  studied  the  dynamic  response  of  bituminous  concrete. 
However,  these  investigators  (U3,  145>  ^6)  have  attempted  to  characterize 
the  viscoelastic  response  of  bituminous  mixtures  by  using  the  phase  angle 
and  the  complex  elastic  modulus,  E  ,  as  obtained  from  the  input-output 
traces  of  dynamic  tests  employing  sinusoidal  loading.  This  investigation 
is  distinguished  by  the  fact  that  a  new  approach  has  been  taken  to  deal 
with  the  sinusoidal  test  results  through  the  use  of  transfer  functions  to 
characterize  the  dynamic  response  of  bituminous  mixtures. 

Form  of  the  Transfer  Function 

From  Tables  3,  k   and  5  it  is  seen  that  the  transfer  function  derived 

from  the  frequency  spectrum  for  each  individual  case  is  of  the  form 

A(s  +  a1)(s  +  a2)(s  +  a  ) 
G(S)  =  (s  +bl)(s  +b2)(s+b3)(s  fbu)         (31) 
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where  A  is  a  constant  and  a.  's  and  b.  's  are  the  corner  frequencies  as 
previously  defined  and  which  are  obtained  in  the  process  of  approximating 
the  frequency  spectrum  with  asymptotes .  The  significant  feature  of  this 
form  of  the  transfer  function  is  that  all  the  factors  in  both  the  numera- 
tor and  denominator  are  to  first  power.  It  may  be  recalled  that  this  is 
due  to  the  geometric  procedure  followed  in  approximating  the  frequency 
spectrum  with  asymptotes  of  slopes  of  six  decibels  per  octave.  It  was 
also  mentioned  that  any  slope  can  be  used,  in  which  case  it  is  obvious 
from  Eq.  29  that  the  transfer  function  may  contain  fraction  powered  or 
higher  powered  terms  in  either  the  numerator  or  the  denominator  or  both. 
This  form  is  not  desirable  since  the  resulting  transfer  function  will  be 
too  difficult  or  even  impossible  to  handle  through  Laplace  transforms. 

From  the  procedure  followed  in  approximating  the  frequency  spectrum 
with  asymptotes  it  is  obvious  that  the  number  of  factors  in  the  numerator 
and  denominator  of  the  transfer  function  in  Eq.  31  depends  on  the  number 
of  asymptotes  used  to  approximate  the  frequency  spectrum.  It  was  found 
from  the  analysis  of  the  experimental  results  that  the  transfer  function 
obtained  using  eight  asymptotes  was  accurate  enough  for  practical  pur- 
poses, which  is  given  in  Eq.  31. 

Check  for  the  Transfer  Function 
From  the  method  of  obtaining  the  transfer  function  discussed  under 
"Determination  of  the  Transfer  Function",  it  is  obvious  that  the  derived 
transfer  function  is  simply  an  analytical  expression  for  the  frequency 
spectrum,  since  the  asymptotic  approximation  is  a  form  of  curve  fitting. 
Hence,  the  absolute  value  of  the  derived  transfer  function  calculated 
at  any  frequency  within  the  spectrum  should  match  the  corresponding  value 
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of  the  gain  or  the  ratio  of  the  magnitude  of  the  output  a.     to  that  of  the 
input  a.,  read  from  the  curve  at  that  frequency,  that  is, 


|G(s)|  =  |G(ja>)|  = 


(32) 


This  has  been  checked  for  all  the  derived  transfer  functions  and  the 
agreement  was  observed  to  be  excellent  for  the  entire  range  of  frequencies 
tested.  A  typical  case  is  illustrated  for  specimen  z-1,  tested  at  80°F. 
Table  11  gives  the  calculated  and  observed  values  which  are  plotted  in 
Fig.  21  as  a  function  of  the  frequency.  Hence,  it  is  of  interest  to  note 
that  the  absolute  value  of  the  transfer  function  at  a  given  frequency  is 
nothing  but  the  dimensionless  output-input  ratio  of  the  system  at  that 
frequency. 

It  was  observed  from  the  experiments  that  at  frequencies  higher  than 
10  rad/sec.  the  output  magnitude  became  too  small  to  measure  (the  higher 
the  frequency  the  smaller  the  output  magnitude).  That  is,  as  the  input 
frequency  cd  approached  infinity  the  output  magnitude  a?  approached  zero, 
from  which  it  follows 


«  0 


(33) 


0)  -*  eo 


where  a  is  the  input  magnitude.  Hence,  from  Eqs .  32  and  33  it  is  at  once 
seen  that 


G(0C3) 


G(s) 


0    (3*0 


0)  -»  co 


JO) 


S  -»  oo 
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TABLE  11 
|G(s) I  and  Output-Input  Ratios 
Specimen:         z-1 
Test  Temperature:  80°F 

r(«\  (0.05)(s  +  0.025)(s  +   Q.2)(s  +  2) 

UKS)   ~      (s  +  0.02)(s  +  0.1)(s  +  l)(s  +  10) 


(JU 


rad/sec 


|G(s) 


a2/al 


0.025 
O.050 


0.0218 
0.0192 


0.0220 

0.0195 


0.100 
0.250 
0.500 


0.0158 
0.0120 
0.0097 


0.0170 
0.0110 
0.0090 


1.000 
2.500 
5.000 


0.0077 
0.0060 
0.0046 


0.0080 
O.OO65 
0.0050 


10.000 


0.0035 


o.ooi+o 


where 


a  =  input  magnitude 
a  =  output  magnitude 
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Eq.  31-*  offers  a  means  of  checking  the  validity  of  the  experimentally  de- 
rived transfer  function,  given  by  Eq.  31 


A(s  +  a^Cs  +  a2)(s  +  a„) 
G(S)  =  (s  +b1)(s  +b2)(s  +b3)(s  +bu)         <31> 


It  is  easily  proved  that  the  limiting  value  of  G(s)  as  s  -»<»  in  Eq.  21 
is  zero. 

Effect  of  Mix  Type  on  Transfer  Function 
The  two  different  bituminous  mixes  studied  in  this  investigation 
were  quite  dissimilar  in  their  aggregate  gradings  (Fig.  10)  and  in  their 
binder  contents,  namely,  5.0£  and  3-5$  by  weight  of  aggregate  for  Mix-1 
and  Mix-2,  respectively.  However,  it  is  apparent  from  Fig.  22  that  there 
was  no  appreciable  difference  in  the  transfer  functions  of  the  specimens, 
irrespective  of  their  composition  or  orientation  when  tested  at  any  one 
temperature.  Even  though  the  grading  and  asphalt  content  were  different 
for  the  two  mixes,  their  density- void  ratio  characteristics  were  similar, 
namely,  1U2.6  pcf  and  G.%   for  Mix-1  and  II42.O  pcf  and  8.5$  for  Mix-2. 

It  can  be  anticipated  that  density  will  have  a  marked  influence  on 
the  response  to  loading  of  a  given  mix  and  hence  on  its  transfer  function. 
However,  more  work  is  necessary  to  study  the  effect  of  this  and  other  mix 
variables. 

Effect  of  Specimen  Orientation 
It  was  pointed  out  under  "Experimental  Investigation",  that  for  each 
mix  three  specimens  were  cut  along  three  mutual  perpendicular  directions, 
designated  x,  y  and  z.  It  is  seen  from  Fig.  22  that  there  is  not  much 
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difference  in  the  transfer  functions  of  the  differently  oriented  speci- 
mens of  either  mix  at  any  one  temperature.  This  suggests  that  effect  of 
anisotropy,  if  any,  is  not  reflected  in  the  transfer  function. 

This  is  contrary  to  the  observations  of  Busching  (35)  who  noted  dif- 
ferences in  response  of  differently  oriented  specimens  prepared  in  the 
same  manner  as  in  the  present  investigation,  though  for  static  loads .  One 
possible  explanation  is  that  the  stress  levels  he  applied  to  the  speci- 
mens, namely,  23-1,  U0.8  and  58.3  psi  were  much  higher  than  the  stress 
level  of  10.0  psi  used  in  the  present  study. 

Effect  of  Temperature 

The  experiments  in  the  present  investigation  were  conducted  at  70°, 
80°,  and  90°F.  By  examining  the  absolute  value  curves  in  Fig.  22,  it  is 
readily  seen  that  temperature  plays  an  important  role  in  the  transfer 
function  of  any  given  specimen.  Asphalt  cement  being  thermoplastic,  the 
asphaltic  concrete  specimen  will  yield  increased  displacements  for  in- 
creases in  temperature  at  a  given  stress  level.  This  results  in  a  de- 
crease in  viscosity  of  the  binder  in  the  asphaltic  concrete  as  tempera- 
ture increases.  Thus  the  output-input  ratio  for  a  constant  input  at  any 
frequency  will  increase  with  decrease  in  viscosity,  or  increase  in  tem- 
perature. This  is  precisely  what  is  shown  by  the  curves. 

Examination  of  the  transfer  functions  in  Tables  3>  ^>  and  5  reveals 
that,  for  a  given  specimen,  the  various  factors  in  the  numerator  and  de- 
nominator do  not  seem  to  vary  much  with  change  in  temperature,  but  that 
the  coefficient  A  varies  appreciably.  It  would  be  of  interest  to  study 
this  change  in  A  with  reference  to  the  viscosity  of  asphalt  cement  or  the 
asphalt  cement-filler  matrix.  This  is  suggested  for  further  study. 
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Effect  of  Specimen  Size  on  Transfer  Function 
Frequency-response  tests  were  conducted  on  three  different-sized 
specimens  with  constant  height-width  ratio  using  Mix-1  at  90°F.  They  were 
all  cut  along  the  same  direction  in  the  compacted  beam.  The  large,  medium 
and  small  sizes  were  1  1/V'  x  1  l/k"   x  2  1/2",  1"  x  1"  x  2"  and  3/k"   x 
3/V'  x  1  1/2",  respectively.  Their  frequency  spectrums  are  shown  in  Fig. 
23,  and  the  transfer  function  is  given  in  Table  12.  The  absolute  value 
curve  for  the  derived  transfer  function  and  the  output-input  ratio  curves 
for  the  three  specimens  are  shown  in  Fig.  2k. 

From  Fig.  23  it  is  seen  that  the  frequency  spectrums  for  the  three 
specimens  are  nearly  the  same  and  that  they  can  be  approximated  by  one 
single  transfer  function.  The  absolute  value  curve  for  this  transfer 
function  and  the  output-input  ratio  curves  for  the  three  specimens  plotted 
in  Fig.  2k   confirm  that  the  transfer  function  is  independent  of  the  speci- 
men size,  since  all  of  them  fall  on  the  same  curve. 

In  the  light  of  this  observation,  it  appears  that  the  transfer  func- 
tion can  be  used  to  represent  the  dynamic  characteristics  of  a  viscoelas- 
tic  material  in  much  the  same  way  as  the  elastic  modulus  represents  the 
stress-strain  characteristics  of  an  elastic  material.   In  other  words, 
for  a  given  viscoelastic  material  such  as  a  bituminous  concrete  there 
will  be  only  one  transfer  function  for  a  given  temperature. 

Prediction  of  Displacements  for  Sinusoidal  Loads 

The  discussions  so  far  have  indicated  the  general  validity  of  the 

transfer  function  for  bituminous  concrete  as  derived  from  dynamic  tests . 

It  was  seen  under  "Related  Theory"  that  once  the  transfer  function  of  a 

system  is  known,  it  can  be  used  to  predict  the  output  of  the  system  for 
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TABLE  12 
Effect  of  Specimen  Size 

LEGEND 
Large:  1  1/V  x  1  l/k"   x  2  1/2" 
Medium:  1"  x  1"  x  2" 
Small:  3/V  x  3/V  x  1  1/2" 
Mixture :  Mix-1 
Test  Temperature:  90°F 

,  x    (0.055)(s  +  0.035)(s  +  0.23)(s  +  2) 
^S;  "   (s  +  0.02)(s  +  0.1)(s  +  l)(s  +  10) 


Calculated  and  Measured  Displacements 
for  Constant  Loads 


Specimen 

f 
o 

P 

X 

(*) 

Calculated 

Measured 

lbs 

psi 

inch(xl0"  ) 

inch(xl0"  ) 

Large 

3-9 

7.8 

15-6 

2.5 

5-0 

10.0 

4.32 

8.6U 

17.28 

4.00 

7.50 

16.25 

Medium 

2.5 

5.0 

10.0 

2.5 
5.0 

10.0 

2.76 
5-52 

11.  ou 

2.25 

3.75 

10.00 

Small 

2.8 

5.6 
11.2 

5.0 
10.0 
20.0 

3.10 
6.20 

12.1*0 

2.00 

3.75 
10.00 
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any  given  input  which  is  a  function  of  time.  After  the  transfer  functions 
for  the  asphaltic  concrete  test  specimens  under  investigation  were  deter- 
mined, attempts  were  made  to  predict  the  output  displacement  of  the  speci- 
mens when  subjected  to  a  sinusoidal  load  input  of  known  magnitude  and  fre- 
quency. Each  specimen  was  tested  under  three  sinusoidal  inputs  at  each  of 

« 

the  three  temperatures  used  previously. 

As  previously  indicated,  the  solution  for  the  output  displacement  of 
a  dynamic  system  with  a  transfer  function  G(s)  for  a  sinusoidal  load  in- 
put is  given  by 

x(t)  =  f  |G(s)|  sin(o3t  +  $)  (20) 

where   f  =  magnitude  of  the  input 

a)  =  frequency  of  the  input 

<t>  =  phase  angle  between  output  and  input. 
In  Eq.  20,  x(t)  is  a  maximum  when  sin(uit  +  4>)  =  1,  so  that 

x(t)    =  f  |G(s)|  (35) 

max    o '     ' 

Experimentally,  it  is  convenient  to  measure  x(t)    at  the  peak  of  the 
sinusoidal  displacement  output.  Thus,  the  calculated  displacement  from 
Eq.  35  can  be  compared  with  the  measured  displacement  for  a  given  magni- 
tude and  frequency. 

It  was  mentioned  elsewhere  that  in  the  experimental  set-up  of  this 
investigation,  the  input  was  recorded  in  pounds  and  millivolts  and  the 
output  was  recorded  in  inches  and  millivolts.  Thus  the  output-input  ratio 
can  be  either  dimensionless  or  in  units  of  inch  per  pound.  The  frequency 
spectrum  and  the  transfer  function  analyses  were  obtained  as  dimensionless 
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quantities.   In  applying  units  to  Eq.  35,  an  experimental  constant  K  needs 
to  be  introduced  and  it  becomes 


*<*>,*«  "Kfo  'G(S)I  (36> 


From  the  sensitivity  control  values  of  the  recorder,  K  was  calculated  and 
found  to  be 

K  =  25  x  10~  inch/pound.        (37) 

The  calculated  and  measured  values  for  the  maximum  displacement  for 
all  the  specimens  are  shown  in  Tables  6,  7  and  8.  In  all  cases  the  mag- 
nitude of  the  sinusoidal  input  was  kept  the  same  and  only  the  frequency 
was  varied.  Comparison  of  the  calculated  and  the  measured  displacement 
values  clearly  indicates  the  close  agreement  between  these  values  in  all 
cases.  Besides  proving  the  efficacy  of  the  transfer  function  as  a  dis- 
placement predicting  tool  for  viscoelastic  materials,  this  also  shows 
that  the  bituminous  concrete  behaves  as  a  linear  system  at  the  levels  of 
stress  considered  in  the  tests. 

Prediction  of  Displacements  for  Static  Loads 
As  was  discussed  under  "Related  Theory",  the  transfer  function  can 
be  used  to  predict  the  time -dependent  displacements  of  a  system  subjected 
to  a  step  function  input  and  a  static  load  can  be  treated  as  a  step  func- 
tion. The  solution  for  the  displacement  was  shown  to  be 


b  t     b  t     b  t     b,t 
x(t)  =  C^e    +  C  e    +  C  e  3     +  C^e    +  C      (27) 


where  C.  's  are  constants  for  the  system  as  defined  earlier  and  b.  's  are 
1  J  1 
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the  roots  of  the  denominator  of  the  transfer  function  which  are  observed 
to  be  negative  in  the  case  of  the  experimentally  derived  transfer  func- 
tions. 

Two  boundary  conditions  are  applicable  to  Eq.  27:  l)  when  t  =  0, 
2)  when  t  =  ». 


When  t  =  0, 


that  is 


When  t  =  oo, 


x(t)  =  0  (38) 


C5  =  Cl  +  C2  +  C3  +  V  (39) 


x(t)  =  C5  (kO) 

When  a  constant  load  is  applied  to  a  viscoelastic  material,  the  dis- 
placement approaches  a  constant  value  after  a  certain  time  depending  upon 
the  nature  of  the  material.   In  the  laboratory,  for  the  bituminous  con- 
crete specimens  tested,  this  time  was  observed  to  be  in  the  order  of  a 
few  minutes.  Hence,  the  second  boundary  condition  discussed  above  can  be 
applied  to  this  steady  displacement  value  in  the  test,  that  is,  Eq.  kO 
can  be  used  to  calculate  the  displacement.  The  limitations  of  using  Eq. 
kO   are  recognized  in  that  in  a  test  t  does  not  reach  »  and  hence  the 
measured  value  will  be  slightly  lower  than  the  calculated  value. 

In  order  to  get  the  calculated  displacement  in  proper  units,  Eq.  kO 
should  be  multiplied  by  the  experimental  constant,  described  in  the  pre- 
vious section  so  that 


x(t)  =  K  C  (Ul) 
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where  K  =  25  x  10   inch/pound .  Cj-  has  the  units  of  pound . 

The  calculated  and  measured  values  of  the  displacement  at  three  tem- 
peratures for  all  the  specimens,  each  under  three  different  constant  loads, 
are  shown  in  Tables  6,  7  and  8.  Comparison  of  the  calculated  and  the 
measured  displacement  values  clearly  indicates  the  close  agreement  be- 
tween the  two  in  almost  all  cases.  This  observation  is  very  significant 
and  it  brings  to  focus  the  following  points : 

1.  The  mathematical  theory  of  transfer  functions  is  applicable  to 
viscoelastic  materials  in  general,  and  to  bituminous  concrete  in  particu- 
lar. 

2.  The  techniques  developed  to  derive  the  transfer  function  in  this 
investigation  are  sound. 

3.  The  transfer  function  serves  as  a  connecting  link  between  the 
responses  of  the  material  tested  under  dynamic  and  static  loads . 

Differential  Equation  from  Transfer  Function 
Once  the  transfer  function  is  known  for  a  given  system,  the  differ- 
ential equation  for  the  system  behavior  can  be  written  using  the  Laplace 
inverse  transform.  For  bituminous  concrete,  the  transfer  function  is  of 
the  form 


A(s  +  a1)(s  +  a2)(s  +  a.^) 
G(S)  =  (s  +bl)(s  +b2)(s  +b3)(s  +bu)         (31) 


From  definition 


G(s)   =     £kl  (io) 

f(s) 
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where  x(s)  =  operational  output  =  <£  (x(t)} 
f(s)  =  operational  input  =  oC   (f(t)) 
From  Eqs.  10  and  31? 

x(s)  A(s  +  a   )(s  +a2)(s+a3) 

?(s)  (s  +  b1)(s  +  b2)(s  +  b3J(s  +  b^)  ^d) 

i 

Rewriting, 

x(s)(s  +  b1)(s  +  b2)(s  +  b3)(s  +  b^)  =  A  f(s)(s  +  8^)(s  +  a2)(s  +  a^    (U2.l) 

Multiplying  and  rearranging  the  terms, 

x(s)[s  +  Bxs3  +  B2s2  +  B  s  +  BjJ  =  A  f(s)[s3  +  B  s2  +  B,s  +  B  ]    (1*3) 

where   B  =  b  +  b  +  b  +  b^  (hk) 

B2  =  bxb2  +  b.b0  +  b3bu  +  bl4b1  +  bub2  +  b3bx         (1*5) 

B3  =  blb2b3  +  b2b3b^  +  b3bUbl  +  Vlb2  {kG) 

Bk  =  \^3\  =    (1.7) 

B  =  aL  +  a2  +  a^  (1*8) 

B6  =  ala2  +  a2a3  +  a3al  ^) 

B?  =  axa2a  (50) 

Rewriting  Eq.  1*3, 

u o p 

s  x(s)  +  B  sJx(s)  +  B  s  x(s)  +  B  sx(s)  +  B.x(s) 

=  A[s3f(s)  +  B 5s2f(s)  +  B6sf(s)  +  B  f(s)]        (^3.1) 
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Applying  the  Laplace  inverse  transform, 


£*l£  ♦  Bx£m  t  BfiM  ♦  B  Mil  ♦  V(t) 


dt '  *     dt3  "     dt" 

"  ft1  t  B^  t  B6-M  .  B/(t)]        (51) 

Eq.  51  is  a  fourth  order  linear  differential  equation  with  constant  coef- 
ficients which  describes  the  dynamics  of  the  bituminous  concrete.  The  in- 
put (force)  and  the  output  (displacement)  are  the  two  variables  in  the 
equation  which  are  functions  of  time.  The  constants  B.  's  are  easily- 
determined  from  the  roots  of  the  denominator  and  the  numerator  of  the 
transfer  function. 

In  dynamical  systems,  if  the  roots  of  the  denominator  of  the  trans- 
fer function  are  real  and  distinct,  then  the  system  is  overdamped  (Hi). 
For  the  bituminous  concrete  under  investigation  the  four  roots,  -b  ,  -b  , 
-b  and  -b,  ,  are  real  and  distinct,  which  indicates  that  the  bituminous 
concrete  is  an  overdamped  system. 

The  constants,  B.  through  B„,   were  calculated  for  each  specimen 
tested  at  each  temperature.  They  are  shown  in  Table  13.   It  is  of  signi- 
ficance to  note  here  that  the  roots  of  the  differential  equation  are  ob- 
tained as  the  corner  frequencies  directly  from  the  frequency  spectrum  for 
the  material  under  test,  and  no  assumption  whatsoever  has  been  made  in 
deriving  the  equation. 
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TABLE  13 
Constants  of  the  Differential  Equation 


Temp 
°F 

Speci- 
men 

A 

Bl 

B2 

B3 

*4 

B5 

B6 

B7 

70 

2-1 

x-1 

y-i 

0.035 
0.024 
0.035 

11.120 

11.114 
12.270 

11.322 
11.266 
23.245 

1.222 
1.166 
5.460 

0.020 
0.015 
0.100 

1.911 
1.800 
3.130 

0.364 
0.323 
1.173 

0.009 
0.006 
0.032 

z-2 
x-2 

y-2 

0.032 
0.035 
0.035 

11.120 
11.120 
11.120 

11.322 
11.322 
11.322 

1.222 
1.222 
1.222 

0.020 
0.020 
0.020 

1.985 
1.785 
1.785 

0.337 
0.300 
O.300 

0.007 
0.006 
0.006 

80 

z-1 

x-1 

y-i 

0.050 
O.OhI 

o.o4o 

11.120 
11.120 
11.120 

11.322 
11.322 
11.322 

1.222 
1.222 
1.222 

0.020 
0.020 
0.020 

2.225 
1.998 

2.485 

0.455 
O.361 
0.489 

0.010 

0.008 
0.009 

z-2 
x-2 

y-2 

0.040 
0.04l 
0.041 

11.120 
11.120 
11.120 

11 . 322 
11.322 
11.322 

1.222 
1.222 
1.222 

0.020 
0.020 
0.020 

2.485 
1.998 
1.998 

0.489 
O.361 
O.361 

0.009 
0.008 
0.008 

90 

z-l 
x-1 

y-i 

0.100 

0.082 

0.100 

11.115 
11.120 
11.115 

11.266 
11.322 
11.266 

1.166 
1.222 
1.666 

0.015 
0.020 
0.015 

1.920 
I.830 
2.190 

0.378 
0.374 
0.383 

0.006 
0.009 

0.007 

z-2 
x-2 

y-2 

0.070 
0.085 
0.085 

11.120 
11.115 
11.115 

11.322 
11.266 
11.266 

1.222 
1.166 
1.166 

0.020 
0.015 
0.015 

2.700 
2.020 
2.020 

0.505 
0.400 
0.400 

0.012 
0.007 
0.007 
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Elastic  Modulus  from  Transfer  Function 
The  differential  equation  that  was  obtained  for  the  bituminous  con- 
crete using  the  transfer  function  relates  the  two  time-dependent  vari- 
ables, force  and  displacement .  If  the  relationship  is  assumed  to  be  time 
independent,  as  in  an  elastic  analysis,  then  the  time  derivatives  of  x 
and  f  drop  off  in  Eq.  51,  and 

BjX  =  AB?f  (52) 

where  x  =  displacement, 

f  =  applied  force, 

A,  Bi  ,  and  B„  are  constants  as  determined  before. 
Eq.  52  gives  the  force-displacement  ratio  as 


I-  k;  <»> 


From  Eq.  53  >  the  stress -strain  ratio  can  be  computed  for  the  material  to 
give  the  elastic  modulus . 

In  this  investigation,  the  gage  length  of  the  specimen  was  two  inches 
and  the  area  of  the  specimen  on  which  force  was  applied  was  one  square 
inch.  Hence,  the  elastic  modulus  E  is  obtained  from 


„    stress     f     .  /,-,,-> 

E  =  itrTin-  *  x72  PS1-  (5I° 


Substituting  Eq.  5**  into  Eq.  53,  and  using  the  experimental  constant  to 
get  the  proper  units, 

E=  KAB7~PSi  (55) 
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where  K  =  25  x  10~  inch/psi.  It  may  be  noted  that  K  is  in  units  of 


inch/psi  in  this  case  because  of  the  area  of  the  test  specimen  being  one 
square  inch  where  the  applied  load  in  pounds  was  considered  as  the  applied 
stress  in  psi. 

It  is  of  interest  to  note  that  Eq.  55  involves  only  the  use  of  con- 
stants  determined  in  the  transfer  function  analysis  for  the  material  under 
test.  Hence,  if  the  differential  equation  derived  from  the  transfer  func- 
tion analysis  is  considered  as  the  viscoelastic  representation  of  the  ma- 
terial, then  it  automatically  includes  the  elastic  case  which  is  only  a 
special  case  of  the  problem. 

Using  Eq.  55  the  elastic  modulus  was  calculated  for  each  specimen 
at  each  temperature.  The  results  are  summarized  in  Table  1*4.  The  values 
are  comparable  to  those  found  in  the  literature  for  bituminous  concrete 
(35>  ^2).  With  respect  to  the  temperature  of  test,  the  elastic  modulus 
decreases  with  increase  of  temperature  which  should  be  the  case  because 
of  the  decrease  of  the  viscosity  of  the  binder  with  increase  in  tempera- 
ture. At  any  given  temperature,  the  orientation  of  the  specimen  does  not 
seem  to  affect  the  elastic  modulus . 

It  was  observed  previously  that  the  mix  type  did  not  seem  to  have 
appreciable  influence  on  the  transfer  function,  due  possibly  to  similarity 
in  density-void  ratio.  The  results  of  Table  1*4  clearly  show  that  this 
holds  for  the  elastic  modulus  also. 


TABLE  Ik 
Values  of  Elastic  Modulus 


Specimen 

E  in  psi 

@  70  "F 

@  bO"F 

@  90 "F 

z-1 

k 
5.05  x  10 

3.20  x  10** 

14 
2.00  x  10 

x-1 

6.68    " 

k.Q2     " 

2.18     " 

y-i 

7.13     " 

k.k5     " 

1.72     " 

z-2 

7.15 

11 

kM 

it 

1.82 

11 

x-2 

7.50 

it 

1+.82 

11 

2.02 

11 

y-2 

7.50 

11 

U.81 

it 

2.02 

it 

Ik 
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CONCLUSIONS 

In  this  study  specimens  of  asphaltic  concrete  cut  from  laboratory- 
compacted  beams  of  two  different  compositions  were  tested  both  statically 
and  dynamically  at  three  test  temperatures.  Based  on  the  results  and 
within  the  limitations  of  this  investigation,  the  following  conclusions 
are  enumerated : 

1.  The  viscoelastic  or  time-dependent  characteristics  of  an  asphal- 
tic concrete  can  be  represented  by  a  transfer  function  G(s)  which  is  a 
function  of  the  complex  variable  s .  The  transfer  function  is  unique  for 
the  material  at  a  given  temperature.   It  is  possible  to  obtain  this  func- 
tion experimentally  from  a  series  of  sinusoidal  load  tests  on  the  material. 

2.  The  transfer  function  for  an  asphaltic  concrete  is  of  the  form 


A(s  +  aL)(s  +  a2)(s  +  a^ 
G(S)  =  (s  +bx)(s  +b2)(s  +b3)(s  +bu)         (31) 


where  A  is  a  constant,  and  a.  's  and  b.  's  are  roots  of  the  numerator  and 

'      l        l 

denominator,  respectively. 

3.  The  roots  of  the  denominator  of  G(s),  namely  b  ,  b?,  b_  and  b,  , 
are  real  and  distinct,  thus  indicating  that  bituminous  concrete  behaves 
as  an  overdamped  system. 

k.     Parameters  obtained  in  the  transfer  function  for  bituminous  con- 
crete are  believed  to  be  better  indicators  of  material  performance  than 
those  commonly  used,  such  as,  Young's  modulus,  Poisson's  ratio,  etc. 


76 


which  change  with  rate  and  time  of  loading. 

5.  Temperature  is  the  one  single  factor  which  has  the  greatest  ef- 
fect on  the  transfer  function  of  asphaltic  concrete.  Increase  in  tem- 
perature increases  the  value  of  the  constant  A  in  the  transfer  function 
equation. 

6.  The  transfer  function  is  a  powerful  tool  useful  in  predicting 
the  displacement  of  asphaltic  concrete  under  an  applied  load,  dynamic  or 
static.  By  treating  the  static  load  as  a  step  function  of  time,  the  re- 
sulting displacement  can  be  calculated  by  means  of  the  transfer  function 
derived  experimental 3 y  from  the  dynamic  test.  The  excellent  agreement  be- 
tween the  calculated  and  measured  values  of  the  displacement  in  this  in- 
vestigation validates  the  concept  that  the  transfer  function  represents 

a  material  property  which  is  independent  of  the  type  of  load  input. 

7.  Through  the  use  of  the  transfer  function  and  without  assuming 
any  spring-dashpot  model,  it  is  possible  to  represent  the  time-dependent 
behavior  of  asphaltic  concrete  by  a  fourth  order  linear  differential  equa- 
tion with  constant  coefficients.  The  coefficients  can  be  computed  from 
the  roots  of  the  denominator  and  numerator  of  the  transfer  function. 

8.  It  is  possible  to  determine  the  elastic  modulus  of  asphaltic 
concrete,  assuming  it  as  an  elastic  system,  through  the  differential  equa- 
tion discussed  above  by  setting  all  the  time  derivatives  to  zero.  The 
range  of  values  thus  calculated  for  all  the  specimens  in  this  investiga- 
tion appears  reasonable  in  the  light  of  values  obtained  by  other  inves- 
tigators . 
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SUGGESTIONS  FOR  FURTHER  RESEARCH 

In  attempting  to  apply  the  concept  of  transfer  functions  to  practical 
problems,  the  limitations  of  this .investigation  should  be  borne  in  mind. 
The  present  investigation  is  a  first  step  in  formulating  and  testing  the 
concept  under  controlled  laboratory  conditions,  assuming  that  the  bitumi- 
nous concrete  behaves  as  a  linear  system.  The  levels  of  applied  stress 
and  the  resulting  small  displacements  were  limiting  factors  in  the  labora- 
tory investigation.  The  conditions  in  full-scale  field  problems  are  con- 
siderably different  from  those  in  the  laboratory  and  it  can  not  be  assumed 
that  the  transfer  function  value  determined  from  limited  laboratory  test- 
ing will  be  specifically  applicable  to  all  field  conditions.  Pilot  inves- 
tigations under  full-scale  field  conditions  are,  therefore,  necessary  to 
establish  the  validity  of  the  transfer  function  concept  for  analysis  of 
problems  of  practical  interest. 

As  a  further  application  of  the  transfer  function  in  the  laboratory 
investigation  attempts  were  made  to  predict  the  displacement  of  the  as- 
phaltic  concrete  test  specimens  subjected  to  shock  loads.  A  weight  of 
500  grams  was  dropped  on  the  specimen  from  a  height  of  four  inches  and 
the  displacement  of  the  specimen  during  the  impact  was  recorded  by  a  Beck- 
man  electronic  recorder  using  an  LVDT  as  the  sensing  device.  Considering 
the  force  due  to  the  impact,  fT,  as  a  Dirac  impulse  function,  the  dis- 
placement was  calculated  from  Eq.  12 
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x(t)  =  c£'{G(s)  f(s))  (12) 

where  G(s)  is  the  transfer  function  of  the  bituminous  concrete  and  f(s) 
is  the  Laplace  transform  of  the  force  due  to  impact,  which  is  f  itself, 
that  is 

f(s)  =  fj  (56) 

Using  the  method  of  partial  fractions  the  soltuion  for  x(t)  was  obtained 
as 


-b  t      -b  t      -b  t      -b,t 
x(t)  =  C^e  x  +  C2e    +  C-e  d  +  Cue  4  (57) 


where  C.  's  are  constants  independent  of  t  and  b.  's  are  the  roots  of  the 
denominator  of  the  transfer  function  G(s). 

Though  the  calculated  and  measured  values  of  the  displacements  agreed 
reasonably  well  in  some  cases,  in  general  the  agreement  was  poor.  This 
was  considered  mainly  due  to  the  lack  of  a  well -designed  shock  loading 
apparatus  capable  of  producing  reproducable  loads.  However,  the  general 
trend  of  the  test  results  strongly  indicated  that  with  further  improve- 
ment in  the  shock  testing  procedure,  better  agreement  might  be  expected 
between  the  measured  values  and  those  calculated  using  transfer  functions. 
This  is  suggested  for  further  research. 

In  this  investigation  attempts  also  were  made  to  measure  the  dis- 
placement of  the  test  specimen  in  the  x  and  y  directions,  both  for  the 
static  and  sinusoidal  load  inputs.  It  was  observed  that  in  a  sinusoidal 
load  test  the  output  displacements  in  the  lateral  directions  were  also 
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sinusoidal,  but  the  magnitude  of  the  output  traces  in  these  directions 
was  too  small  for  meaningful  interpretation.   However,  by  the  use  of  an 
electronic  recorder  with  greater  amplification  than  the  one  used  in  this 
study  it  should  be  possible  to  get  enlarged  output  traces.  Then  from  the 
input-output  data  of  the  test  specimen  a  transfer  function  could  be  de- 
rived for  each  direction,  say  G(s)  ,  G(s)  and  G(s)  ,  the  input  being  the 

x      y        z 

same  for  all  the  three  cases .   It  would  be  of  interest  to  study  the  re- 
lationship among  these  three  transfer  functions  in  order  to  throw  more 
light  on  the  use  of  transfer  functions  to  characterize  material  behavior. 
This,  again,  is  suggested  for  further  research. 

In  addition,  the  following  suggestions  are  offered  for  extension  of 
this  investigation  to  studies  which  will  contribute  understanding  for  the 
use  of  transfer  functions  in  the  structural  design  of  pavements. 

1.  The  effects  on  the  transfer  function  of  viscosity  of  the  asphalt 
cement  and  the  density  of  a  bituminous  mix  need  to  be  studied  in  detail, 
since  this  would  lead  to  a  better  understanding  of  the  viscoelastic  char- 
acteristics of  the  mix  itself,  and  thus  possibly  to  a  rational  mix  design. 

2.  The  concept  of  transfer  functions  can  be  extended  to  a  layered 
system  such  as  an  asphaltic  concrete  beam  or  a  slab  supported  on  springs 
or  on  actual  base  course  materials.  The  laboratory  evaluation  of  the 
transfer  function  for  such  a  system  would  pave  the  way  for  a  full-scale 
pavement  analysis  in  the  field. 

3.  Soil  and  concrete  are  the  other  pavement  materials  to  which  the 
transfer  function  analysis  can  be  logically  extended.  Since  the  Theo- 
logical behavior  of  clay  soils  has  been  studied  widely,  it  would  be  of 
interest  to  apply  the  transfer  function  concept  to  them.   It  is  probable 


80 


that  their  mechanistic  behavior  could  be  better  understood  through  the 
use  of  transfer  functions. 

k .  The  analysis  in  this  investigation  has  been  based  on  the  assump- 
tion on  linearity.  However,  pavement  materials  do  not  always  conform  to 
this  assumption,  and  hence,  it  would  be  logical  to  study  the  nonlinear 

i 

transfer  functions  in  these  cases.   Such  an  analysis  should  place  the 
structural  design  of  pavements  on  a  more  rational  basis. 
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APPENDIX 
TYPICAL  EXAMPLE  OF  OBTAINING  THE  TRANSFER  FUNCTION 

Typical  test  results  in  a  sinusoidal  loading  test  on  an  asphaltic 
concrete  specimen  (z-l)  at  80°F  and  the  frequency  response  curves  plotted 
from  them  are  given  in  Table  10  and  Fig.  20,  respectively. 

The  first  step  is  to  approximate  the  frequency  spectrum  with  asymp- 
totes of  slopes  of  six  decibels  per  octave  and  to  note  the  corner  fre- 
quencies, as  shown  in  Fig.  25  and  Table  15. 

TABLE  15 
Slopes  of  Asymptotes  and  Corner  Frequencies 


Slope 
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The  approximate  modified  transfer  function  was  seen  to  be 

1                        2     1                              n     n-1 
—r  — -r —  £- 

G(jco)   =  A(jcu  +  o^)        (jco  +  o>2)  ...    (jco  +  con)  (28) 

Substituting  the  values  from  Table  15  in  Eqn.   26,  we  get 

_,.    n  A(jco  +  0.025)(ja>  +  0.2) (joo  +  2.0)  ,cQv 

G(jCD)  "     (Jcd  1  0.02)(Jcd  +  O.lHjco  ;  1.0)(ja)  ;  10.0)  <58) 

and  the  transfer  function  is  accordingly 

C(s)    ■        Ms   +  0.025)(s   +  0  2)(s   +  2.0)  .      . 

^S;   "      (s   +  0.02)(s   +  0.1)(s   +  1.0)(s   +  10.0)  ^y; 

where  the  constant  A  is  determined  at  one  experimental  frequency  as  shown 
below . 

A  frequency  in  the  spectrum  is  chosen  at  which  the  experimental  curve 
and  a  horizontal  asymptote  intersect  -  such  as  to  =  O.k   rad/sec  in  the 
present  case.  From  Eqn.  58,  the  absolute  value  of  the  transfer  function 
is  obtained,  by  virtue  of  the  property  of  a  complex  number,  as 

|G(d»)|  =     A(cu2  +  0.  OOO625  )1/,2(cd2  +  O.OU)1/2^  +  l4.0)l/2 

(co2  +  0.000U)1/2(cd2  +  0.0l)1/2(o)2  +  l)1/2(co2  +  100)1/2 

(60) 
Substituting  for  to  =  0.U  in  Eqn.  60,  and  evaluating  we  obtain 

JG(oai)  J^  =  Q>l4  =  0.20U  A  (6l) 

From  Fig.  25  the  value  of  the  gain  at  o>  =  0.^  rad/sec  is  read  as 
-  40.0  decibels;  that  is  experimentally 
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°^2^  =  gain  =  -  Uo.O  db  (62) 

input 


Converting  Eqn.  6l  to  decibel  units  and  equating  it  to  Eqn.  62, 

20  log10(0.2C4A)  =  -  UO.O  (63) 

Solving  for  A  in  Eqn.  63,  A  is  obtained  as  a  dimensionless  constant, 

A  =  0.05  (Sk) 

Substituting  for  A  in  Eqn.  595  the  transfer  function  is  obtained,  as 

nf    \    (0.05)(s  +  0.023)(s  +  0.2)(s  +  2.0)       ,£- 
G(s)  "  (s  +0.02)(s  to:i)(»n!o)(8  f  10.0)    (65) 
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